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An oscillator with stochastic frequency is discussed as a 
model for evaluating the quantum coherence properties of a 
physical system. It is found that the choice of jump statistics 
has to be considered with care if unphysical consequences are 
to be avoided. We investigate one such model, evaluate the 
damping it causes, the decoherence rate and the correlations 
it results in and the properties of the state for asymptotically 
long times. Also the choice of initial state is discussed and its 
effect on the time evolution of the correlations 

PACS numbers: 03.75.-b, 42.50.Lc, 03.65.-w, 05.30.Ch 



I. INTRODUCTION 

The harmonic oscillator plays a pivotal role in the de- 
velopment of physics, both classical and quantum me- 
chanical. It can be solved exactly, it displays a variety of 
non-trivial features and it has been applied successfully 
to a variety of phenomena. 

A straightforward generalization of the simple oscilla- 
tor is one with time dependent frequency; its Hamilto- 
nian is 



2m 2 w 



(1) 



With a deterministic time evolution of u>(t), this can de- 
scribe parametric frequency generation, radio-frequency 
traps or frequency modulated signal transmission. If, on 
the other hand, the frequency u)(t) is behaving in a ran- 
dom way, the model can be applied to a variety of phys- 
ical situations, vide infra. The oscillator with a random 
frequency, has been the subject of a variety of theoretical 
approaches .see e.g. the review [Q. In quantum theory, 
the model (□) can be taken as an approximation to a sit- 
uation where the oscillator is coupled to some degrees of 
freedom carrying out their quantum dynamics in a com- 
plicated fashion. Then it may be possible to model this 
by a time evolution so irregular that it can be regarded 
as stochastic. 

With all quantum processes, we are not only inter- 
ested in the time evolution of the averaged variables and 
their noise characterized by the dispersion from the av- 
erages. We also want to know how coherences are con- 
served, which variables retain their quantum character 



and how the classical behavior emerges. In this paper we 
undertake the task to investigate these questions for the 
oscillator, when oj(t) is given by a stochastic process. 

Stochastic processes can be of two major types: In the 
one ui(t) varies continuously according to some model 
of the Brownian motion type or a modification of this. 
The ensuing theory is mathematically well understood, 
amenable to analytic approximations and often physi- 
cally transparent. In another class of models, random 
telegraph signals called, the variable ui(t) jumps between 
constant values according to some stochastic rule. In sim- 
ple cases, the consequences of such models can be eval- 
uated 0, they have simple properties and they can be 
used as approximations for continuous behavior in the 
proper limit [^|^]. In a quantum treatment, they offer 
the advantage that between the jumps, the well known 
properties of the harmonic oscillator can be applied di- 
rectly. Thus we choose to utilize such a model here, and 
evaluate the influence the stochastic evolution has on the 
quantum coherence properties. 



II. PHYSICAL BACKGROUND 

One application of a stochastic oscillator is to the 
s.c. continuous Stern-Gerlach measurement discussed by 
Dehmelt and his colleagues |j [7j]. Here a trapped parti- 
cle, an electron, sees two different trap potentials depend- 
ing on the internal spin state. This is changing due to its 
own dynamics, and by monitoring the oscillational fre- 
quency of the charge, the experimentalist can see the fre- 
quency follow the changes of the internal spin state. This 
constitutes a continuous observation of the two quantum 
states, and the coherence between them cannot survive 
for any observable time, and a well defined frequency 
is always observed. This continuous Stern-Gerlach ef- 
fect should be described by a stochastic quantum model 
along the lines developed by Carmichael and colleagues 
||[)| . They have been able to demonstrate how the mon- 
itoring of the state of a two level system gives precisely 
the type of well resolved jumping which is characteristic 
of the continuous Stern-Gerlach behavior. We have at- 
tempted such modelling, but the computational capacity 
needed was found to be unreasonably large, and hence we 
have settled on the stochastic process described above. 
This should suffice when one wants to investigate the 
consequences for the survival of quantum coherence in 
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the oscillator. 

Another situation where the time evolution is inter- 
rupted by change of the potential function governing the 
motion is given by molecular dynamics on adiabatic en- 
ergy surfaces. Here spontaneous emission transfers a 
wave packet from an excited state to the ground state 
[ fTpf or laser coupling may transfer the state back and 
forth between the potential levels. If this process is co- 
herent, Rabi flopping may be observed jll]] or adiabatic 
motion may ensue |12| . With less coherent coupling, the 
dynamic evolution has been described by surface hopping 
Jl3| , where the deterministic motion on one potential sur- 
face is stochastically transferred to another one. It is well 
justified to describe this process by a random switch be- 
tween the various potential curves. In this paper, the 
model assumes that the two harmonic states are situated 
directly above each other; this may not be a realistic de- 
scription of a molecule, but it simplifies the treatment, 
and does not greatly affect the features we are investigat- 
ing. It would not, in principle, be difficult to carry out 
the simulations on more realistic modifications of this 
model. 

When an evolving quantum state is suddenly trans- 
ferred to another potential at random times, we expect 
the quantum coherence to be obliterated, and the sys- 
tem to behave in a more classical way. This is seen in 
the density matrix in such a way that it tends towards the 
diagonal both in the position and the momentum repre- 
sentation. The uncertainties along the diagonal are then 
large, their product is expected to greatly exceed the min- 
imum allowed by quantum theory. We have performed 
calculations to observe such disappearence of coherence. 

The simplest model is to let the stochastic frequency 
jump between two values u)i and LO2 with the constant 
jump probability v. This model will be found to have 
physically less satisfactory features: The coherences de- 
cay, but the jumping process feeds energy into the system 
so efficiently that the energy grows at a rate quite com- 
parable to the rate of vanishing of the coherences. The 
explanation for this strange behavior is in the physics of 
the oscillator motion. The system spends much of its 
time near the classical turning points of the oscillation, 
and consequently the jumps most likely take place here. 
Whenever the state jumps from the turning point of the 
flatter (slower) potential to the steeper (faster) one, it 
appears far up the potential slope, and a large increase 
in energy ensues. Jumping in the opposite direction, the 
system looses energy, but not enough to compensate for 
the gain. Thus the rapid increase in energy. 

To overcome the unphysical features of our model, we 
strive to conserve energy in the jump process. This means 
that we have to favor jumping in the region where the 
change of potential energy is minimal, i.e. near the bot- 
tom of the wells. Thus the wave function should be most 
eager to jump, when it overlaps most with the ground 
state; then only the kinetic energy is transferred between 
the levels. This also maximizes the conservation of lin- 
ear momentum; near the potential minimum, the motion 



most resembles free propagation. 

In order to achieve the goal outlined above, we choose 
the following model: When the state on level 1 is ^i(t), 
we jump to level 2 with the probability 



Pia(t) = i/|<*§ |*!(t)> 



(2) 



where ^\ i s the ground state of level 2. The jumps are 
taken to occur at a random time determined by Pn, and 
after that the state *£i(t) is transferred to level 2 and 
evolved according to the corresponding potential. The 
probability to jump back to level 1 is now determined 
by P2i(t) given by the expression (0) with 1 and 2 in- 
terchanged. This simulation continues, and enough nu- 
merical data are accumulated to allow us to evaluate the 
ensemble averaged quantities. The results are then com- 
pared with those obtained by the simple model having 
P = v constant. 



III. TIME EVOLUTION OF THE STATE 

The Schrodinger time evolution with the potential ([!]) 
is given by 

d . s h 2 d 2 V(x,t) 1 2 2 



ih— *(x,t) = - 



2m + 5 ™^*(M), (3) 



where i = 1 or 2. Using the standard definition of the 
Wigner function JlJ] 

! f x xf Px 

W{X,P) = —jd X (X + -\p\X-~)e^(-i- 

(4) 



we can obtain the equation of motion 



dW{X,P) ( P\ dW(X,P) 



at 



771 J 



dX 



dW(X, P) 
dP 



= 0. 
(5) 



For the harmonic oscillator, this can be solved, and the 
various moments evaluated. It is, however, easy to see 
that the equations of motion for the moments close in 
each order, and their evolution equations can be obtained 
by simple calculations from that of the Wigner function. 
The mixed moments of the type (xp) denote symmetrized 
products, as these are given by the Wigner function. 

The equations for the first moments we find to be the 
classical dynamic equations 








For the second order ones we find 




(6) 
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Given the initial values, it is straightforward to integrate 
these equations for the moments. From these we can 
calculate the variances 
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(p 2 (t)) - (p(t)) 2 

(x(t)p(t)) - (x(t))(p(t)). 



(8) 
(9) 

(10) 



If the initial quantum state is a Gaussian wave function, 
it suffices to integrate the moments to the first two orders; 
it is known that the state will stay Gaussian, and that 
it is described by its moments of the first two orders. A 
straightforward calculation shows that the Schrodinger 
equation is satisfied by the state 



^>{x,t) = {2^l{t))- 1 / 4 eM- 



{x-(x(t))f 
4*2(t) 



2ha 2 (t)' 



{x{t))f 



*<P(*)> 



(x-(x(t)))+t6(t)), (11) 



where the phase 9{t) disappears in the calculation of the 
Wigner function (0) and thus does not affect the physics. 
All moments occurring in the state are not independent. 
If we start from a minimum uncertainty state, the com- 
bination 



(12) 



is conserved. 

The density matrix in the position representation is 
given by (x\ \ p \ x 2 ), and its quantum coherence is 
the amount of off-diagonality, i.e. the dependence on 
the variable x = x\ — x 2 . Its moments can be evaluated 
directly from the Wigner function according to 



dx\dx 2 {xi - x 2 ) n (xi I p \ x 2 ) 



2Tr(ih) n J dR 



gn 

~dP 



^W(R,P) 



(13) 



p=o 



A similar result can be obtained for the off-diagonality 
in the p-representation. 



IV. NUMERICAL RESULTS 

Starting from a Gaussian initial state, the equations 
(|J) and §j§ are integrated, and the state can be directly 
obtained from Eq. (O) . We have found this integration to 



be orders of magnitude faster than the direct integration 
of Schrodinger's equation, and obtaining the state from 
( |Tl| ) , we can directly calculate the Wigner function. The 
most time consuming part is still to accumulate enough 
single histories to obtain the ensemble averaged Wigner 
function and the corresponding average moments. This 
computation is, however, straightforward, and consists in 
simulating histories with the proper jump probabilities 
given as explained above. The whole calculation is time 
consuming but essentially trivial. 

In the numerical work we want to evaluate the average 
position and momentum, and their variances. We have 
used dimensionless scaled units, which bring the numeri- 
cal results into ranges of order unity. We have integrated 
the above problem using a Gaussian initial state on level 
1 with the initial values 



(x(0)) = 2.0 
a 2 x (0) = 0.5 



(P(0)> = ; 
a 2 (0) = 0.5 ; <(0) = 0. 



(14) 



The stochastic model consists of jumping between the 
two frequencies lui = 0.7 and lo 2 — 1-2 with the jump 
frequency parameter v = 0.8. We use two models, in the 
one the jump probabilities are constant 



P 12 = P 2 i = 



(15) 



in the other one we use the nonconstant, time dependent 
jump probabilities determined as explained in Sec. 2 and 
Eq.(||). In order to get satisfactory statistics for eval- 
uating the ensemble averages, the number of individual 
histories making up the ensemble consists of N — 30000 
runs. 

In Fig. 1 we report the energy of the oscillator as 
function of dimensionless time. The solid line is the con- 
stant jump probability result, which displays a nearly 
exponential growth. The dotted line is the result with 
the nonconstant jump probability determined according 
to Eq.(|]). We can see that the goal declared has been 
achieved; the energy grows much more slowly and clearly 
not exponentially. This is ascribed to the tendency of the 
nonconstant jump model to concentrate the jumping to 
the center of the potential. 

That this is, indeed, the case can be seen from Figs. 2 
and 3. In Fig. 2 we show the statistical distribution of 
jumps from the level 1 to 2. The classical turning points 
of the oscillational motion are indicated by the solid bars 
on the horizontal axis. The solid line shows the result for 
the constant probability, which leads to many jumps near 
or outside the turning points, where the oscillator spends 
large times. The peak at the right hand edge indicates 
the position of the starting wave packet. The dashed line 
is the result of the simulation with nonconstant probabil- 
ities. The jumps occur much closer to the center of the 
potentials as expected, even if the effect of the turning 
points is still seen. The total jump rate has decreased, 
thus the initial asymmetry is no longer so manifest. In 
Fig. 3, the distribution of jumps from state 2 to state 1 is 
shown. The constant probability case shown by the solid 
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Ensemble energy 




FIG. 1. The ensemble energy as a function of time with 
constant and nonconstant jump probability. The jump fre- 
quency is v — 0.8. The initial wavefunction is Gaussian 
with parameters (X(0)) = 2.0, <P(0)) = 0.0, (cr|(0)) = 0.5, 
(cr|,(0)> = 0.5, (iTxp(O)) = 0.0. The frequencies are o>i = 0.7 
and u>2 = 1-2. The initial frequency is u\. The ensemble used 
is of size N = 30000. 



line still concentrate to outside the turning points in the 
steeper potential; these are again shown as solid bars. 
The symmetry is now in the opposite direction because, 
with the initial condition chosen, there must occur one 
initial jump to state 2 before we can have a jump back. 
The dashed line shows the result for the model with non- 
constant probabilities. The same features are seen as in 
the previous figure, but in this case of faster oscillations, 
the turning points have slightly enhanced effect. The de- 



Number of jumps 



constant 
nonconstant 




FIG. 3. The number of jumps from the potential ui2 to wi. 
The parameters are the same as in FIG. 1 



Number of jumps 




FIG. 2. The number of jumps from the potential with u>\ to 
the potential with oj^. The bars on the X-axis show the clas- 
sical turning points with the frequency wi. The parameters 
are the same as in FIG. 1 



crease of energy growth in the case of nonconstant jump 
probabilities allows the model to damp the average mo- 
tion faster. This is seen in both the position variable, Fig. 
4, and the momentum variable, Fig. 5. These variables 
are expected to behave in very similar manner for the 
harmonic motion. The damping time can be estimated 
to be of the order of 10 units, which is much less than the 
one observed with constant probability. For harmonic 
motion, the second moments of both position and mo- 
mentum are expected to grow like the energy, see Fig. 1. 
This is also seen in the variances; for constant probabil- 
ities they increase exponentially, while the nonconstant 
case grows more slowly. This is shown in Fig. 6, which 
should be compared with Fig. 1; the momentum vari- 
ance behaves in a very similar manner. For the model 
with nonconstant jump probabilities, the time scale of 
motional damping has been found to be of order 10, the 
range over which energy grows is several times this value, 
see Fig. 1. If we look at the dephasing causing the disap- 
pearance of quantum correlations in the position variable, 
we plot the expectation value 

{x 2 ) — J J dxidxi (xi - x 2 ) 2 (xi | p | x%), (16) 

which can be obtained directly from the Wigner function 
as shown in Eq. (|l3|). The position density matrix can 
also be obtained from the Wigner function by inverting 
the relation (|J) 
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Variance of position 



Ensemble position 



constant 
nonconstant 




FIG. 4. Ensemble position as a function of time. The pa- 
rameters are the same as in FIG. 1 



Ensemble momentum 




FIG. 5. Ensemble momentum as a function of time. The 
parameters are the same as in FIG. 1 




FIG. 6. Position variance as a function of time. The pa- 
rameters are the same as in FIG. 1 



(xi | p | x 2 ) — j dP exp 



iP(x 1 - x 2 ) 



W p 



(17) 



Fig. 7 shows the decay of the quantum mechanical corre- 
lations, and we can see that they disappear after a time 
of approximately 4 units, which is considerably faster 
than the other time scales. In particular, the increase in 
energy is totally negligible over this time for the noncon- 
stant jump probability; see Fig. 1. This clearly corrects 
the unphysical behavior of the model with constant prob- 
ability. The constant model, on the other hand, destroys 
the quantum coherence about as fast as the other model 
does, but the asymptotic value for large times does not 
seem to vanish but linger on at a small but finite value. 
However, the very process of jumping removes the pos- 
sibility to retain coherence, the actual statistics of the 
process seems to matter less. 

When the harmonic oscillator is subjected to a random 
perturbation, we expect the final result to resemble a 
thermal distribution. In order to check this, we plot the 
density matrix in the occupation number representation. 
This is directly obtained from the position representation 
( |l7| ) using the relation 



{tii | P | n 2 ) 



dx 1 dx 2 {ni | xi)(xi \ p \ x 2 )(x 2 \ n 2 ). 

(18) 



The states (x | n) are the Hermit polynomial eigenstates 
of the harmonic oscillator, and the integrations in Eq.jT^) 
are thus straightforward. 

The diagonal elements of the density matrix in the oc- 
cupation number representation with the frequency u>i 
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function in the X, P— plane. In Fig. 9, the original Gaus- 
sian state at t = 0, rotates and distorts due to the influ- 
ence of the jumping, which here is taken with noncon- 
stant probabilities. The fore front rotates faster and the 
distribution loses phase information within a few units 
of the time. In Fig. 10, the time variable is 30, and the 
phase information has disappeared. This picture corre- 
sponds to the same situation as that in Fig. 8, and the 
dip seen in the middle corresponds to the minimum at 
n = in that figure. It is obvious, that, in the occupation 
number representation, there are no off-diagonal density 
matrix elements in the state shown in Fig. 10. 



Wigner function at t^O.O 



Wigner function at t=3.0 



FIG. 7. X-coherence as a function of time. The parameters 
are the same as in FIG. 1 



are shown in Fig. 8 for both models at the asymptoti- 
cally large time 30. We find that the model with noncon- 
stant probability decreases much faster and goes to zero 
for large values of the quantum number. Both follow 
asymptotically an exponential decrease with the quan- 
tum number n. 



Wigner function at t=6.0 



Wigner function at t=9.0 



Diagonal elements 




FIG. 8. Diagonal elements of the ensemble energy density 
matrix (m\g\n) as a function of time. The Fock-states have 
the frequency wi . The scale on the y-axis is logarithmic. The 
parameters are the same as in FIG. 1 

The exponential decrease indicates that the state is 
close to a thermal distribution of the Planck type. How- 
ever, the states near n = clearly deviate more than 
the following ones from a thermal distribution. This is 
verified when we follow the time evolution of the Wigner 



FIG. 9. The Wigner function at four different times with a 
nonconstant jump probability. The parameters are the same 
as in FIG. 1 



Wigner function at t=30.0 




FIG. 10. The Wigner function at t — 30.0 with a noncon- 
stant jump probability. The parameters are the same as in 
FIG. 1 
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V. EFFECT OF INITIAL SQUEEZING 



The numerical results obtained so far have started from 
the ideal minimum uncertainty state (0). In order to 
see how sensitive our results are to this assumption, we 
have integrated the equations also with squeezed initial 
states. In the X-squeezed state we have set cr^-(O) = 0.25 
and t7p(0) = 1.0 ; in the P-squeezed state we have 
0-2.(0) = 1.0 and crp(0) = 0.25. Fig. 11 reports the 
growth of energy in the model with nonconstant jump 
probabilities; compare this figure with Fig. 1. For the P- 
squeezed state, the advantages of letting the jump proba- 
bility vary is partly lost; the X-squeezed case is not very 
different from that shown in Fig. 1. The wave packet 
still spends much time near the turning points, and in 
the case of P-squeezing at the initial point, the wave 
packet is very broad here. This increases the overlap 
with the central region of the oscillators, and the jumps 
adding energy to the system become more numerous. In 



Ensemble energy 



Number of jumps 





X squeezed 

— P squeezed 







FIG. 11. The ensemble energy as a function of time 
with X-squeezed o"x(0) = 0.25, crp(0) = 1.0 and P-squeezed 
CT x(0) = 1-0, o"p(0) = 0.25 as initial conditions. The other 
parameters are the same as in FIG 1. 

Fig. 12 we show the jump distribution from the initial 
level 1 to level 2, and we see that the P-squeezed state 
clearly tends to concentrate around the turning points. 
This figure should be compared with Fig. 2, and we can 
see that even the X-squeezed state has an enhanced jump 
probability as compared with the unsqeezed result, the 
dashed line in Fig. 2. Fig. 12, does not, however, show 
any asymmetry as found in the model with a constant 
probability. The jump frequency from level 2 back to 
level 1 shows a behavior almost identical to that of Fig. 
12; the P-squeezed state tends to jump around the turn- 
ing points, the distribution has a minimum at the center, 
which is not found in the case of X-squeezing. Fig. 13 
shows how the average position is decaying. There is lit— 



X Squeezed 
P squeezed 




FIG. 12. The number of jumps from the potential with wi 
to the potential with u)2- The bars on the x-axis show the clas- 
sical turning points with the frequency u)i. The parameters 
are the same as in FIG. 12 



tie difference between the P- and X-squeezed cases, and 
only a slight decrease in the decay rate is seen when com- 
pared with the decay in Fig. 4. The average momentum 
behaves in a very similar way. 



Ensemble position 




FIG. 13. Ensemble position as a function of time, 
parameters are the same as in FIG. 12 



The 



The variance in the position is shown in Fig. 14. For 
the case of X-squeezing, the wave packet is broad near 
the center, where the jumps take place, and the jumping 
consequently adds little width to the packet; the behavior 
is very similar to that of the dashed curve in Fig. 6 show- 
ing the result for the unsqueezed case. The P-squeezed 
case, however, jumps more often near the turning points, 
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and consequently the width in position space is growing 
more rapidly. However, not by far as rapidly as with 
the constant jump probability shown in Fig. 6, the solid 
curve. When we look at the disappearance of the quan- 



Variance of position 



X Squeezed 
P squeezed 




FIG. 14. Position variance as a function of time. The pa- 
rameters are the same as in FIG. 12 





X Squeezed 

— P squeezed 









FIG. 15. X-coherence as a function of time. The parame- 
ters are the same as in FIG. 12 



turn coherence, the variable fllq), the result in Fig. 15 
shows that squeezing does not affect the rate of disap- 
pearance of quantum coherence much as compared with 
Fig. 7. Jumping always destroys the coherence rapidly 
and efficiently as expected. 



VI. CONCLUSION 

In this paper we have discussed the time evolution of 
quantum coherence and correlations in a harmonic os- 
cillator with a stochastic frequency. For simplicity we 
choose a jump process, where the frequency stays con- 
stant but switches randomly according to a stochastic 
rule. In Sec. II of the paper we discuss the possible physi- 
cal situations which may be modelled by such behaviour. 

We find that one should be careful in modelling ran- 
dom jumping, because the enforced transitions violate 
energy conservation, which leads to clearly unphysical 
growth of the oscillator energy. In order to amend this 
shortcoming, we device an ad hoc model, where jumping 
near the potential minimum is enhanced, which tends to 
conserve energy. In addition, the motion resembles free 
propagation most closely at this point, and hence the 
jump achieves momentum conservation as close as pos- 
sible, when we jump here. This gives some additional 
physical motivation to our model, when we consider it 
as an attempt to emulate the behaviour of real physical 
processes. 

The model suggested is found to achieve the asserted 
goal: the energy grows but little, and the decoherence 
time scale is much shorter than that over which the en- 
ergy changes. We find that jumping in itself tends to 
destroy quantum coherence, but jumping near the cen- 
ter of the potentials achieves this without a build-up of 
unphysical energy. After an asymptotically rather long 
time, the diagonal elements of the density matrix resem- 
bles that of a thermal state, but the occurrence of very 
low quantum numbers is less likely than it would be ther- 
mally. This is verified by plotting the Wigner function, 
which displays a ring shaped form resembling that of an 
operating laser. No phase information survives for large 
times. 

We have always started the integration near one of the 
turning points. Squeezing the state in the P-direction, we 
find that the accompanying broadening in the X-direction 
leads to more jumps near the turning points. This trends 
to counteract the improvements introduced by our non- 
constant jump probability model. With squeezing in the 
X-direction less ill effects are found. In this case, they 
would, however, be seen if we started the wave packet 
closer to the center of the potial. 

We conclude that the model of a stochastic oscilla- 
tor can be used to investigate the quantum coherences, 
but care has to be excersised when the model is chosen. 
In our work, we have shown how to achieve physically 
reasonable results with a suitably modified jump rate. 
Which model is appropriate in a given physical system 
must be left open; different situations may require differ- 
ent choices of jump statistics. 
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